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Abstract 

The 'Miranda Procedure' proposed for analyzing Dalitz plots for CP asymmetries in 
charged B and D decays in a model-independent manner is extended and refined. The 
complexity of CKM CP phenomenology through order A 6 is needed in searches for New 
Dynamics (ND). Detailed analyses of three-body final states offer great advantages: (i) 
They give us more powerful tools for deciding whether an observed CP asymmetry rep- 
resents the manifestation of ND and its features, (ii) Many advantages can already be 
obtained by the 'Miranda Procedure' without construction of a detailed Dalitz plot de- 
scription, (iii) One studies CP asymmetries independent of production asymmetries. We 
illustrate the power of a second generation Miranda Procedure with examples with time 
integrated rates for B^jB^ decays to final states Kstt + tt~ as trial runs with comments on 
B ± K ± 7i + 7v-/K ± K + K~. 
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1 Landscape of B Uj d jS , D u ^ s & r CP Violations 

The predictions of CKM theory have been impressively confirmed to a degree that has 
persuaded a part of our community to focus on scenarios of Minimal Flavour Violation 
(MFV) - i.e., models of New Dynamics (ND) that contain the same sources of flavour 
violations as the Standard Model (SM). Some intriguing work has been done along these 
lines, yet we view the hypothesis of MFV as far from compelling at present. The data 
still allow for sizable deviations from SM predictions in heavy flavour transitions; in D° 
and some B s decays they could even be dominant. Furthermore, baryogenesis requires 
the intervention of ND with CP violation. If ND appears around the (9(1 TeV) scale 
underlying the weak-electric phase transition, which is not intrinsically connected with 
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flavour dynamics, it will affect CP asymmetries in heavy flavour decays - but not on the 
leading level. CP asymmetries are given by three classes of observables, namely 

• \q/p\ 1j which shows purely indirect CP violation. 

• Absolute amplitudes \Af\ ^ \Af\ that show purely direct CP violation and depend 
on the final states. 

• The relative phases between q/p and Af <g> A%, which will depend on the final state 
/. We write Af <g> A*j rather than just AfA*f, because for a three-body final state 
one has to denote the position in the two-dimension plot. The significance of this 
feature will become clearer through the illustrations given below. 

1.1 Present Status of CP Asymmetries in B, D Sc r Transitions 

Oscillations have been observed for all three heavy flavours B s , B d and D° mesons, but 
on very different numerical levels. 

1.1.1 Indirect CP Violation in B ds Transitions 

Indirect CP violation has been measured with very good accuracy in B d — > i{}Ks/Kl and 
B d -> 7T + 7T~ [I]: 

S(B d ^ipK s ) = 0.658 ±0.024 (1) 
S(B d ->■ 7T+7T-) = -0.61 ±0.08 (2) 

Purely indirect CP violation gives S(B d — > ipKs) = —S(Bd — > 7r + 7r _ ). 
Very recent data from LHCb [2] on B s — > ipf Q {980) give 

<p s = -0.002 ±0.083 ±0.027 rad (3) 

1.1.2 Direct CP Violation in B U)d)S Transitions 

Direct CP asymmetries has been established in B d decays [HE]: 

A CP (B d ^ K + n-)\ PDG , w = -0.098 ±0.013 (4) 
A C p(B d -> K + 7r")| L HCb'ii = -0.088 ±0.011 ±0.008 (5) 
C{B d -> TT+TT - ) = -0.38 ±0.17 (6) 

No sign for direct CP violation has been found in B d — > Ksir - by sizable asymmetry 
can still be allowed: 

A CP (B d K S 7T°) | PDG ao = 0.00 ± 0.13 (7) 
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Intriguing evidences for direct CP violation has been found in B s and B + in (quasi-)two- 
body final states [H EJ E] : 

A C p(5 s ^# + 0| L HCb'ii = +0.27 ±0.08 ±0.02 (8) 

A CP {B + -»■ D CPl+1] K + ) = +0.24 ±0.06 (9) 

A CP (B + -> p°K + ) = ±0.37 ±0.10 (10) 

A CP (B + ->■ / (1370)tt + ) = ±0.72 ±0.22 (11) 

-»■ 77^+) = -0.37 ±0.09 (12) 

Acp(£+ -> / 2 (1270)if + ) = -0.68+°;}? (13) 

A CP {B S # + 0|LHCb'ii = ±0.27 ±0.08 ±0.02 (14) 

Acp(B s ->• if + 7T)lcDF = ±0.39 ±0.15 ±0.08 (15) 

1.1.3 Evidence for CP Asymmetries in D° Decays 

No sign of indirect CP violation has been found in D° — > K + K~ /7r + 7r~ [5]: 



0-88i^ , 4>= (-10.2^)° (16) 



Direct CP asymmetry has been found in T(D° — > K + K ) — T(D° — > 7t + tt ) with 3.5 
sigma away from zero by LHCb and with 2.7 sigma by CDF; 

AA CP = -0.82±0.21(stat) ±0.11(syst)% LHCbg] (17) 
AA CP = -0.62±0.21(stat) ±0.10(syst)% CDF[7] (18) 

with AA C p = A C p(D° -»■ K + K~) - A C p(D° -»■ tt+tt"). This is the first significant 
evidence for CP violation in AC ^ dynamics, and it is important whether it is due to 
alone SM or need impact from ND. 

1.1.4 Evidence for CP Asymmetries in r Decays 

In t~ — > Kg-n^v decays one has a prediction [8 J 

A C p(t + + K s ir + )\sM = (0.36 ± 0.01)% , (19) 

independent of dynamics that generate K° — > K° oscillations, and data from the BaBar 
Collab. 0: 

A CP (r + ->■ V ± K S 7l + {> 07T ])| B aBar = ("0.36 ± 0.23 ± 0.11)% . (20) 

1.2 CKM Matrix Parametrization through 0(X 6 ) 

The usually applied the Wolfenstein parameterization of the CKM matrix gives real parts 
through O (A 3 ) and the imaginary part through 0(X 4 ). The CKM matrix is usually 
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described by the four parameters A, p, rj and A with the last three ones of order unity; 
thus one gets iKfc/Kbl ~ Aa/p 2 + rf ~ 0(A). 

PDG states \V ub /V cb \ ~ 0.085. The global fit leads to p ~ 0.13 and 77 ~ 0.34. It means 
that one has to use a parametrization with through order of A 6 and with other quantities 
of true order of unity. One has been found in Ref.[10j: 

/l_Ai_Al_A^ A h\ 4 e- i5 ^ \ 

2 £J 16 ' ' ' 

-A+f/ 2 , 1 - f - ^(1 + 4/ 2 ) - fhX 5 e~ iS ^ fX 2 + hX 3 e~ iSQM 

+ g(4/ 2 -4/i 2 -l), -f/wr^, +C(A 7 )(21) 

/A 3 , -/A 2 - /iA 3 e- l<5 Q M 1 - - fhX 5 e~ iS ^ 

v +*/+*/. -** a ; 

A global fit of the CKM matrix gives: A ~ 0.225, / ~ 0.75, /i ~ 1.35 and the 'maximal' 
phase 5qm — 90°. 

This pattern is not so obvious as from the Wolfenstein parametrization, more subtle 
for CP violation and is similar only in a semi-quantitive way [UJ. To give three examples: 

• CP asymmetry in B d — > ipKs depends in SM on 

_ V3Y* VcbV^ ysin^ QM +(y) 2 sin2^ M 

v tb v t * d v; b v cs 1 + (ff + ?f cos6QM 

One gets: 

S(B d -> tpK s ) = sin20! ~ 0.63 - 0.69 for S QM - 75° - 90° (23) 

S(B d ipK s ) = sin20i ~ 0.74 for 5 QM ^ 100° - 120° ; (24) 

i.e., CKM dynamics could produce S(B d — > ipKs) ~ 0.74 as largest value for CP 
asymmetry with 5qu — 100° — 120°, not with the maximal 5qu = 90°. 



Im 



2(i//ii'w.ti[i//HM 

l + (4ft//)Acosi Q M V ' 



Again one finds that indirect CP violation in B s is CKM suppressed in the SM by 

' y t tVtsVc b v: s 

v tb v t * s v; b v cs 

with S QM ^ (75 - 120)°. 
• Direct CP violation in B — > D + K ± depends on sin0 3 , where one gets: 

, ( VubVud \ n \2/n^ X sin^ Q M , , 

Thus 

03 = 0.28 / 0.34 / 0.42 for 5 QM = 75° / 90° / 110° . (27) 

Therefore sin20i ~ 0.69±0.06 and sin03 ~ 0.34±0.07 are consistent with CKM dynamics 
with lower values of 0i & 03 correlated with each other with 10 % vs. 20%. Nevertheless 
the impact of ND can 'hide' in predicted CP asymmetries. 
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1.3 Present Resume on CP Asymmetries in Two-Body Final 
States in B and D Decays 

SM generates indirect and direct CP asymmetries in B & D (and in K) transitions. Their 
strengths are based on several items: 



• The CKM matrix is discussed above in Sect.lJZ We can say that SM is at least 
the leading source of CP violation in B on most transitions - except at present for 
B s — > ^)0/^/o(98O). On the other hand ND can affect CP asymmetries on the level 
of ~ 10 — 20% for B u d decays. Furthermore one has to focus on correlations with 
CKM suppressed decays of B u ^,s (and K and D( s )) on the level of 20 %. Therefore 
one need more accuracy from data and their interpretation to find impact of one 
(or two) ND - and to probe three-body final states. 

• While the final states K + n~ in B^ and B s are the same, the underlying dynamics 
are very different: 

— SM amplitudes for B^ — > K~tt + are given by 'tree' Cabibbo suppressed tran- 
sitions b — > uus and (1-loop) 'Penguin' b — > sqq; ~ 10 % CP asymmetry seems 
a reasonable value in SM. 

— On the other hand SM amplitude for B s — » K + tt~ is given by 'tree' Cabibbo 
favoured b — > uud and Cabibbo suppressed 'Penguin' b — > dqq. Therefore one 
expects CP asymmetry on the 'natural' level of (9(1%). 

ND can enhance 'Penguin' amplitude significantly. However one expects such 
effect in other B decays - unless a nearby resonance can affect mostly B s , but 
not Bd,u decays. 

For D decays the landscape is much more subtle, but also very 'topical': 

— CKM dynamics produce direct CP asymmetries in singly Cabibbo suppressed 
(SCS) decays around the scale of 0.001. 

— CP asymmetries in double Cabibbo suppressed (DCS) ones are zero at (9(A 4 ). 

• CP asymmetries are controlled by non-perturbative QCD. 

• The difference between Aq P {t + — > v + -Ks7r + [> 7r°])| measure d and Acp(t + — )• v + 
Ks^ + )\sm depend on our control of non-perturbative QCD - like the plots of the 
final states of A CP (r + + [Kn] + ) and A cp (t + [Kim} + ) [TJ]. 

Finding impact of ND in CP asymmetries in B and D decays is one thing - however 
probing the 'shape' of one (or two) ND is another challenge. 

Dedicated studies of three-bodies final states are needed to identify important features 
of ND involved [33] . Three-body final states analyses are very time consuming. In Sect(2] 
we list the general advantages that such analyses merit the needed work; we first sketch 
in Sect|3]the situations for three-body decays of B s , B^ D° in general for searching CP 
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asymmetries through (partly) time integrated data to set the stage; afterwords we give 
more realistic situations in Bd transitions in Sect, ^2 and comments on B s transitions in 
Sectj5] finally we summarize our main conclusions in Sectj6j 



2 Advantages of Studies of Three-Body Final States 

While the weak dynamics from CKM and ND are the driving forces for CP asymmetries, 
one has to control FSI from non-perturbative QCD not only in a qualitatively way - that 
is the focus of our study 

2.1 Opportunities Offered by Dalitz Plot Studies 

No study of any three-body final states in B decays have found an established CP vi- 
olation, and none from K or D mesons shows any sign for it. However there can be - 
actually they are more likely to be found. The average over CP asymmetries in a Dalitz 
plot is expected to be much larger than 'local' asymmetries, which often compensate with 
each other. As explained in some detail in [13] , crucial insights into CP odd dynamics will 
be learnt from their impact on final state distributions. Dalitz studies will play a central 
role in the future for several reasons: 

• Differential or 'local' asymmetries could be considerably larger than ones averaged 
over the Dalitz plot. 

• For two-body final states there is only one CP asymmetry, namely T(B°/D° — > 
h + h~) vs. T(B°/D° — > h + h~). On the other hand the topologies of Dalitz plots for 
B° — > h + h~h° and B° — > h + h~h° are in general different; for example the two half 
of the plots Sh+h° ~ s h~h° f° r the B° and B° are different. However their sum have 
to be symmetric - unless CP asymmetries occur! 

• While the difference in total rates for B/D — > 3h vs. B/D — > 3h are affected by 
production asymmetries, differences between corresponding regions in the Dalitz 
plots are not. 

• Nontrivial correlations provide powerful validation tools. 

• The pattern of a CP asymmetry that has emerged in a Dalitz plot can tell us about 
the spin structure of the underlying effective operator. 

• The cleanest experimental sign whether an observed asymmetry is produced by 
direct or indirect CP violation (or which parts are due to one or the other) is their 
dependence on the time of decay. Direct asymmetry is independent of the time of 
decay, whereas indirect violation evolutes in time in a clear prescription, since it is 
driven by oscillations. With only time integrated data with two-body final states 
one cannot decide it. If one observes a CP asymmetry in a leading CKM final state 
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- like Bd — > ipKg or B s — > ip(f) - you will argue that it is most likely indirect. Yet 
for CKM-suppressed decays you hardly have such an argument. As explained later, 
one can use more decision criterions from three-body final states. 

• Time depending CP asymmetries give us more informations about underlying dy- 
namics for B° and D° transitions. Of course one needs more statistics. Partially 
time integrated rates for three-body final states give us more insights. 

Dalitz studies offer also a more technical advantage when searching for CP asymme- 
tries, namely 'tunable' strong phases. Since that is of direct relevance for our subsequent 
analysis, we explain it next. 



2.2 Phases with Breit-Wigner Resonances 

With CP violation being expressed by a complex weak phase due to CPT invariance, 
it can lead to an observable asymmetry only if one has the interference between two 
different amplitudes. Yet more is needed: hadronization has to affect the two amplitudes 
differently. This is usually expressed by stating that the two amplitudes have to exhibit 
different weak as well as strong phases. 

Consider the decay P — > f receiving contributions from two coherent amplitudes, 

A(P^f) = e^ e^ sl \AA + e^ e^ S1 \A 2 \ (28) 
A(P^f) = e~ l < e^lA.l+e-^ e i5 ™\A 2 \, (29) 

where (j>Y and 5f sl are the weak and strong phases, respectively, and Ai are the moduli 
of the amplitudes. The CP asymmetry between partial widths 

_ T(P -> /) - T(P -> /) 

Acr ~ T(P^ S ) + T(P^!Y (30) 

is given by 



2sin(A0 


w ) sin 


(A5 fsi )|^a: 




1 + \A 2 Ai\ 


2 + 2\A 2 A 1 


cos(A0^) cos(A5 FSI ) 



A CP = 1 • \ A A 19 , ol A A I ~ — ( A XW\ « ( A JCFST\ ( 31 ) 



CP violation is induced by A(p w , but it becomes observable only if the final state inter- 
action (FSI) introduces a non-trivial phase shift. 

For two-body final states it is often implied that the strong phases Sf SI carry a fixed 
value for a given final state / (the two amplitudes are assumed to have different isospin 
contents up to isospin violation). 

In the case of three-body decays, the transition P — > f is dominated by resonant 
intermediate states. The requirement of non-trivial strong phase different is satisfied by 
the energy dependent phases of the resonances. The Breit-Wigner excitation curve for a 
resonance R reads 

Fr W (s) = — - „ , (32) 

m\-s- vm R T R (s) 
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introducing a sizable phase as expressed through 



R K ) (m\- s) 2 + (m R T R (s)y 1 ; 

where r R (s) denotes the energy dependent relativistic width. 

For P P1P2P3 we define = (pi + P2) 2 and s 2 = (pi + P3) 2 . The previously 



constant strong phases and amplitude moduli in Eq.(31) now depend on the position in 
the Dalitz plot, S FSI —> 5(si,S2) and A — > A(si,S2). The resonant amplitudes populate 
the whole phase space in D decays, and a large portion of it in B decays. Therefore, the 
CP asymmetry will also depend on the Dalitz plot coordinates, Aqp — > Acp(si, S2). 

After taking the modulus square of these amplitudes one reads off that a CP asym- 
metry will arise, when there are non-zero weak phases. Having to deal with non-uniform 
strong phases might appear as a complication that just creates more work for analysis of 
decays with three-(and four-)body final states. However there is an award for extra works: 
a resonance - in particular if it is relatively narrow - can tell us where CP asymmetries 
have to be and that the asymmetry has to change over a relatively narrow range in the 
Dalitz plot. The Dalitz plots carry the same area independent of production asymmetries; 
yet the relative corresponding population densities probe CP invariance. 



2.3 'Miranda Procedure' 

Indirect CP violation in B° decays is and will be well measured in B^ — » ipKg and 
B s — > ifxp, ip fo(980) - and maybe in D° — > <pKs\ ND can impact those transitions in a 
sizable way. Direct CP violations affect final states in different strengths, in particular 
CKM suppressed ones both in SM and ND. The existence of ND might - and probably 
will - be found in indirect and direct CP violation in two-body final states - yet its main 
features have to extracted from three-(and four-)body final states. Some asymmetries can 
tell us about the spin operators creating about etc. There ere three classes of sources of 
CP asymmetries, namely 

1. from CP conjugated quasi-two-body final states; 

2. interference between quasi-two-body final states; 

3. contributions from "true" three-body final states or broad resonances like a. 

Contributions from the first class like Kp or ixp are obvious. However from the second 
class one finds positive and negative contributions to the CP asymmetry; therefore those 
get washed out from the total integral over the phase space. This applies mostly to the 
third class of CP asymmetries. Therefore we will denote CP asymmetries from the first 
and second classes by CPVa and CPVi below. For the third class one can help sizably 
from future theoretical efforts. 

The advantages listed above for CP studies in three-body final states justify the con- 
siderable 'overhead' in statistics and tuning efforts that constructing a satisfactory Dalitz 
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model requires. Yet one cannot count on obtaining a unique Dalitz model even with 
infinite statistics. It is thus of considerable practical value to develop another method 
for analyzing a Dalitz plot that is model-independent; it will allow important statements 
about the existence of a CP asymmetry and its approximate localization inside the Dalitz 
plot with smaller data sets than constructing a full fledged Dalitz model. At the very 
least it would help to identify the sub-domain in the Dalitz plot where one had to focus 
the tuning efforts for the Dalitz model. 

In Ref. [13] we have proposed one method that can serve such a purpose in a quantify- 
ing way. When searching for CP asymmetries we have suggested to analyze the significance 

which amounts to a standard deviation for a Poissonian distribution rather than the 
customary fractional asymmetry 

= N(i) - N(i) 

in particle vs. anti-particle populations N(i) and N(i) for each bin i, respectively. One an- 
alyzes, whether the S(i) distribution has a higher frequency of exhibiting large deviations 
from zero than expected for fluctuations We have illustrated this method - 'mirandizing' 
in vernacular - by applying it to B ± — > K ± 7i + tt^ and D ± — > 7r ± 7r + 7T _ ; the Dalitz plots 
have been constructed with fast MC simulations making specific assumptions about the 
underlying dynamics and the source of the CP asymmetry. In those pilot studies we 
could show that using the observable E(i) instead of A(z) allowed a much more robust 
extraction of the seeded CP asymmetry and its location inside the Dalitz plot. It remains 
to be seen of course how mirandizing holds up when treating real data. 

Our claim is not to replace full fledged Dalitz plot analysis. Our goal is to present 
an analysis that can produce significant results on CP violations from smaller data set 
while maintaining many of the advantages of a full Dalitz plot study. The latter is still 
the final goal of our road to 'Rome' - the impact of New Physics on CP violations in 
nonleptonic decays of beauty and charm hadrons. We also want to encourage others to 
try other possible roads to this 'Rome'. An interesting work can be found in Ref. p2 



2.4 Comment on CP Asymmetries in r — > v[Ktt / Kmr] 

The SM generates a global CP asymmetries in r~ — > uKs[S = 0] due to K° — K° 
oscillations with 2Re ejc, but not beyond that. On the other hand ND has a larger chance 
to appear in the CP asymmetries in r~ — > v[Kix'sY , since SM amplitudes are Cabibbo 
suppressed. 

For r~ —7- i>[K7r]~ one has a three-body final states in general, and in t~ — > v[K2tt/3K\ 
one has three-body hadronic final states. The 'Miranda Procedure' can and should be 
applied here with some refining Dalitz plots: the masses of the hadronic systems are not 
fixed as in B and D decays. 
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3 Setting the Stage for Probing CP Invariance 



Both indirect and directly CP violations have been established in K° and Bd transitions, 
and SM through CKM theory gives at least the leading contributions. The ND expected 
- or hoped for - to find around the few TeV scales should produce some 'footprints' 
through CP asymmetries in B and D transitions. However 'reading' them produces 
large both experimental and theoretical challenges. The time evolutions of the flavour 
tagged transitions provide the most powerful tool in identifying to source of the underlying 
dynamics. To begin this projects we want to show what you learn from non-flavour tagging 
transitions without time resolved analyses. 

For this study we describe the time integrated rates also for non-flavour tagged data 
and give comments on partly time resolved rates. 

CP asymmetries are control by five observables: 

1. x — AM/f and y = Ar/2f , which are insensitive to CP violation; 

2- \q/p\ 7^ 1) which shows purely indirect CP violation, and is determined by \q/p\ — 
1 — |a§£, where denote the CP asymmetry in semi-leptonic decays in 'wrong'-sign 
leptons; 

3. absolute amplitudes \Af\ ^ \Af\ that show purely direct CP violation and depend 
on the final states and 

4. the relative phases between q/p and Af <g> A* f) which will depend on the final state / 
due to direct CP violation. We write Af® A* f rather than just A f A*j, because for a 
three-body final state one has to denote the position in the two-dimension plot. The 
significance of this feature will become clearer through the our illustrations later. 

We give general expressions with these observables. Then we show that in describing for 
Bd decays we can ignore yd effects, while for D° decays one has to include both xd and 
yo dependences, but only to first order. For B s transitions one has x s ^> y s , but one has 
to include y s effects due to spectacularly fast x s oscillations for time integrated data; CP 
asymmetries controlled by Im^Aj <g> Af are suppressed by l/x s . 

Indirect CP violation affects all channels through two quantities, namely \q/p\ and the 
relative phase between q/p and Af®Ay, their weight of course depends on \\q/p\ — 1| and 
the strength of As mentioned above we can use the approximation of \q/p\ = 1 

for Bd and B s channels; the effect of indirect CP violation is affected the direct CP on 
rate Bd and B s modes and therefore the impact of ND that is probably different for mode 
to mode. 

For two-body final states like B s — > h + h~ vs. B s — > h + h~ time dependent CP 
asymmetries are reduced by l/x 2 s ; however for direct CP violation in B s — > h + h~h° 
vs. B s — > h + h~h° one can find an asymmetry between corresponding regions of the sum 
of Dalitz plots of B s — >■ / and B s — >■ / due to interference effects - i.e., without flavour 
tagging. 
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Obviously flavour tagged time resolved analyses bring the largest information about 
the underlying dynamics; our main goal for this study is how many lessons can be obtained 
from non-flavour tagged time integrated data. Flavour-tagged and time resolved data will 
come later. 



3.1 Three-Body Decays for Neutral Mesons 



Or 7T + 7T 



it 



For neutral B or D decays into two-body final states like K + K 

that it is a (even) CP eigenstate. For / = h + h~h° like K s ii + 7i~ , K S K + K or n ' n 7T 
the judgement is more complex: it can be [CP=+] K s f (980)/K s a — > K s ir + n^, [CP—-] 



is clear 
o 



A 



or [CP= 



0-7T U ->■ 



7T + 7r 



~7T°. At the same 



K S (j) ->■ K S K + K~, [CP=-j p"?r — y n * n tx 

time one has final states K* ±r K T , p ± 7r =F etc. and interferences between them. Not only 
the total time integrated widths for P — >■ h + h~h° vs. P — > h + h~h° give us a lesson on 
CP violating dynamics, but also their 'topologies' - i.e., the distributions over the Dalitz 
plots. Therefore we denote Af and Aj for P — > h + h~h° and P — > h + h~h°, respectively. 

The 'Miranda procedure' can be applied to all three-body final states, but here we 
will discuss it only for / = h + h~h° like Ks^ + ^ ■ 

The time dependent rates can distinguish the three types of CP violations: \q/p\ ^ 
\p/q\, \Af \ 7^ \Af \ and Im^Aj® AJ ^ Im^Af <8> Ay. For practical reasons we focus on time 
integrated widths for this study: 
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CP violation can appear in the time integrated widths in principle from the three sources 
listed above. 

While the time integrated observables have the largest statistics, the time dependent 
ones have most detailed information about the underlying dynamics; partially integrated 
one can give us most of that information depending on how the parameters for the hadrons 
and their transitions. 



dtr(p^f-t) 

rr = J dtF(P -+ f;t) ; 
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therefore 
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3.1.1 B d -»■ h + hrh° 



For B d transitions one can simplify these expressions in two ways, namely y d — and 
\qd/pd\ — 1 for realistic experimental sensitivities. Integrated over all times of decays we 



get 



dt\A(B d ->f;t)\ 2 -. 

3 _ _ 
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For equal productions of B d and -B^ - p{B d ) 
1 



d 

p{B d ) 
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we get 
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Im 
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+ ; 9 lm 



dt\A(B d ^ f;t)\ 2 + / dt|^(B d ^/;t)| 



l^/| 2 + l^/l 



(44) 



For a two-body final state / = h + h~ one gets no information about CP violation from 
the time integrated sum as expected. However for / = h + h~h° direct CP violation can 
produce an asymmetry in corresponding regions of the sum of Dalitz plots as sketched 
before. Our studies given below will illustrate this feature. Their sum is weighted by 
their ratio due to a real production asymmetry (or a difference in their efficiencies) 



dt\A{B d ^ f;t)\ 2 + (l-p(B d )) 



1 + x 2 



dt\A(B d ^ f;t)\ 2 
'id 



(2p(B d ) + xj)\A f \ 2 + (2(1 - p(B d )) + x 2 d )\A f \ 2 + 2x d (l - 2p(B d ))lm —Af ® A 
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(45) 



If there is production asymmetry, it is not a 'vice', but a 'virtue'. It can be tracked 
by B d — > ipK~7[ + vs. B d — > ipK + 'K~ or by B — > ipK ± . The strength of indirect CP 
violation is measured in B d —> 4>Ks, whether the SM produces the whole or just the 
leading source of it. Used as an input for B d — > K s it + -n~ , K S K + K~ one can interpret 
the impact of direct CP violation through the term Im^Aj® A*^j, see Eq.(45). 

3.1.2 B s -»■ h + h-h° 

For B s transitions, one can assume \q s /p s \ — 1, since even with sizable ND contributions 
to B s — B s \q s /p s can differ from unity not more than several permil. While AT S is small 
- y s — Ar s /2r s ~ 0.094 ± 0.024 - it should not been ignored. Integrated over all times 
of decays we get 
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Therefore we get for non-flavour tagged sum 



roo roo _ _ 

/ dt\A(B s -± f;t)\ 2 + / dt\A(B s ^f;t)\ 
.Jo Jo 
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If there is a production asymmetry one gets: 



roo roo _ _ 

/ dt\A(B s ^f;t)\ 2 + / dt|.A(S a ->/;t)| 

JO JO 



= 2|^| 2 + 2|A / f + 2y s ((\A f \ 2 + l^f) - 2Re (^A f ® A*X) + 2(1 2p(jBs)) Im ® A 



While the strength of indirect CP violation has not measured in time-resolved data on 
B s — > ^0, there are some evidence that it might be significantly larger than the CKM 
prediction of around sin2/3 s ~ 0.03 — 0.05. That situation should be more clarified in 
one to three years. The 2y s term can give us useful information about the dynamics of B 
decays, but itself does not represent a CP asymmetry. In principal if there is a production 
asymmetry - it could be tracked by the Cabibbo suppressed transition B s — > ■ifjK + n~ vs. 
B s — > ipK~n + - one could obtain A f ^> A*^j ; however it is greatly suppressed in 

B s — > K s tt + it~, K S K + K~ by l/x s - i.e., the spectacularly fast oscillation. 



3.1.3 Comments on B° transitions 

Our goal for B° transitions is to analyses the impact of ND on direct CP asymmetries in 
three-body decays in CKM suppressed channels. 

Indirect CP violation affects all transitions of a given meson - Bj, B s and D° - in 
the same way. For Bd we have measured it with good accuracy in Bd — > ipKs with CKM 
dynamics as the leading source. The SM prediction tell us that \qd/Pd\ can differ from 
unity by less than 0.001. 

For B s transitions some evidence has been found in B s — > ipifi and B s — > l^X processes 
for a large impact of ND. We expect that evidence will be validated or reject with good 
accuracy in the foreseeable future from LHCb, CMS and ATLAS. For the time being one 
can use two scenarios, namely 

• Case CKM: sin2/3 s ~ 0.03 - 0.05, \\q s /p s \ - 1| < 0.0001; 

• Case CKM + ND: sin2/3 s ~ 0.11 ±0.02, \\q 8 /p,\ - 1| ~ 0.003 

keeping in mind that such cases will be decided about future data on B s — > ip<p and 
B s ->• l~DX. We consider the impact of ND in B s ->• K S K + K~, K s n + n-. 
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3.1.4 D° -> h+h-h° 

For D° transitions both xd and yn are small and probably of similar size, but \qd/pd\ 
could differ from unity by up to 30 %. Integrating over all times t of D° decays one gets 

roc poo _ _ 

/ dt\A(D° ->• f;t)\ 2 + dt\A(D° ^ f;t)\ 2 cx 
Jo Jo 

\A f \ 2 + \A f f - y D Re (j + A } ® A) + * D Im (j _ ^ A} ® A} (50) 

An production asymmetry leads to 

/ dt\A(D°^f-t)\ 2 + (l-p(D )) dt\A(D° f;t)\ 2 oc 

JO JO 

p( J D°)|A / | 2 + (1 - p(L>°))|A 7 | 2 - y D Re ^p{IT)- + (1 - p(^ u ))- | A T ® + 

+Xi)Im 



(51) 



In D° decays the interplay of indirect and direct CP violations is not so clear mostly due to 
very slow oscillation. Therefore we will consider scenarios with \\qo/pD\ — 1| — 0.1,0.03, 
|t4/|/|Aj| ~ 0.1, 0.03 and the relative phase of Qd/pd and Aj®A* f . We will present studies 
in a future paper. 

3.1.5 Comments on ND scenarios 

For analyzing CP asymmetries in h + h~h° final states one has to include not only PV 
final states, but also SP final states like scalar a and k resonances. One reason for that 
is that exchanges from charged Higgs fields will introduces CP asymmetries already on 
the zero-loop processes, and they will affect scalar resonances more than PV final states. 

3.2 Future Progresses in Describing Dalitz Plots 

When LHCb and the B factories at SLAC and KEK was planned and approved CKM 
theory had a competition with other models for the leading source of CP violation in 
heavy flavour transitions. BaBar and Belle have found with great success that CKM 
provides at least the leading source of the establish CP asymmetries in Bd decays. LHCb 
is in a great position to find whether CKM is also the leading source of CP asymmetries 
in B s decays even for B s — > i/xf) transitions that are largely reduced in CKM. 

In addition LHCb and Super-Flavour Factories have to deal with the difficult task to 
find non-leading source(s) of CP asymmetries in suppressed decays. There are several 
candidates for that ND - even there is no 'standard' version of SUSY, let alone for other 
NDs. 

The 'Miranda Procedure' can allow us to find a clean evidence for a CP asymmetry 
without a theoretical input. It will encourage much more theoretical progress on our 
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understanding on soft QCD, which can make use of other theory tools obtained from 
hadronic dynamics. The tasks one faces in B and D decays into three-body final states are 
not quite as challenging as for the astronomers mentioned above: we know the locations 
where clear CP asymmetries can occur in B° and D° transitions - in particular in Ks^ + ^ 
and K$K + K~ final states: They get sizable contributions from K$p°, K^tt^ and Ks4>- it 
has been shown that the Breit-Wigner parameterization provides a good approximation 
for vector mesons like p, K* and <fi. However scalar resonances will in general not be 
described that way; for a final state Ksf o(d80) it might give a decent description due 
to its relatively narrow width, but not for i^5<r(600) or kit due to their wide widths. 
Still using a Breit-Wigner parametrization published data include true scalar resonances 
under 'non- resonance' label; more theoretical analyses including the treatment of chiral 
dynamics is needed, very topical - and possible now based on progress in the last few 
years. We know that they have to be performed separately for different B/D — > 3h 
transitions. We should understand the following: if present data for K7r + 7i~ final states 
are best fitted without any a — > 7r + 7r~ contribution or for two different cr 12 — > 7i + tt~ one, 
one should not ignore one with the usual single <r(600) as long as it gives a satisfactory 
description. One needs more experimental and theoretical analysis. 

There are several important reasons to analyze the production of scalar resonances in 
the detailed way. Let us just sketch one: Many models of ND contain physical charged 
Higgs states that can introduce CP asymmetries even through their tree-level exchanges. 
Obviously scalar Higgs exchanges will leave their 'footprint' in the production of scalar 
resonances with more weight than for pseudoscalar and vector states. Therefore such ND 
will produce more 'readable' impacts in CP asymmetries with scalar resonances and their 
interferences with pseudoscalar-vector final states. Therefore we can first focus on the 
known location of the peaks of the p, K* and the /o(980) and their widths from available 
data. Using this general input from theory we can generate binning for B° — > Ksh + h~ 
and do it separately for Bd and B s transitions. 

4 Second Generation 'Miranda Procedure' 

The procedure given in Ref. [T3] is obviously powerful for finding CP asymmetries and 
even 'localizing' them. However one wants to make it more quantitatively and to un- 
derstand its source(s) - in particular for B transitions CKM gives sizable 'backgrounds' 
when searching for ND. 

The goal is to find a way to evaluate the strength of local effects, namely to have 
numbers that are equivalent to the asymmetry between (time) integrated rates, Ac?- 
The key idea is to divide the combined Dalitz plot into bins with equal populations. If 
one knows the number of bins where CP is violated, the one can compute a local average 
value of t4cp(si, S2), since the number of events is proportional to the number of bins. 

Each bin has N = N + + N~ events, with N + and N~ being the numbers of B and B 
candidates. We assume that there are regions in the Dalitz plot with at least a few tens 
of bins in which positive events (iV + ) occur with the same probability p. N + follows a 
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binomial distribution with expected value and variance given by 



and 



E[N + ] = Np, V[N + ] =Np(l- p) (52) 

When N is large enough (at least a few tens of events), the Central Limit Theorem 
ensures that N + follows a normal distribution, allowing one to write exact expressions for 
the expected difference N + — N~ . In this case one has 

N+ _ n~ 2N + 

App = — = — - 1 (53) 

CP N N y ' 

with 

f I = E[^] = 2 ^P-l=2p-l (54) 

^ - VMS = «^ = Wp> (8B) 

If CP symmetry is conserved - hereafter we assume that there is no other source of 
charge asymmetry - the probabilities of positive and negative events are equal, p = 1/2. 
One therefore has 

A* = 0, a 2 = 1 (56) 

When CP is violated the Dalitz plot will have regions with and without asymmetries. 
Therefore, the distribution of A^p will be a superposition of a Gaussian with /i — and 
a = 1/ plus some other function representing the CP violating bins. The form of the 
latter depends on how CP violation occurs in the Dalitz plot and also on the specific final 
state. 

Three-body final states result, in general, of a cascade process in which the heavy 
meson decays to a resonance plus a 'bachelor' hadron. The decay amplitude of a heavy 
meson P is usually modeled by a coherent sum of resonant amplitudes, weighted by 
constant complex coefficients 

M = Y^CiAi, Ci = CLie 5 \ (57) 

i 

CP violation results in a difference between JA(P) and M. = A4(P). More specifically, 
M and M may 

• differ by the magnitude of a set of resonant modes 

• or a difference between their relative phases 

• or a combination of both. 

Rescattering at the hadronic level is a long distance effect that mixes different final states, 
e.g. KKix — t> 7r7T7T, and this is another source of CP violation. 



17 



In this paper, we consider three possibilities: (i) CP violation due to re-scattering as 
a constant excess of one specie over the other limited to some region of the Dalitz plot; 
(ii) CP violation through a difference in the magnitude of a resonant amplitude; (iii) a 
difference between relative phases. 

The case of constant CP violation is the simplest: it leads to asymmetries that have 
always the same sign. An uniform and localized excess of one charge state over the other 
is equivalent to a constant value of p. In this case the distribution of CP violating bins 



will be a Gaussian with mean and sigma given by Eqs.(54) and (55). Integration over the 
phase space results is an observable global A CP . 

Differences in magnitude would also correspond, in principle, to a constant p. How- 
ever, the net effect depends on the resonance spin and on the contribution of the other 
resonances. Angular momentum conservation constrains the angular distribution of the 
decay products. For vector particles, for example, the Breit-Wigner is modulated by a 
spin amplitude which is proportional to the cosine of an helicity angle. In the region where 
the momentum configuration is such that the helicity angle is 90°, the amplitude goes to 
zero. The relative contribution of the CP violating amplitude varies from bin to bin, in 
spite of the constant difference in its magnitude. Even in the case of a scalar resonance 
(constant spin amplitude), one needs to take into account the contribution from other res- 
onant amplitudes to the CP violating bins, which in general is not constant. Also in this 
case an integration over the phase space results in an observable global A C p, although 
the local effect will be always diluted by the relative contribution of the CP violating 
amplitudes. 

Differences in phases are the most complex case. As illustrated in [13], such differences 
lead to asymmetries that change sign across the Dalitz plot. Integration over the phase 
space could result in a null asymmetry, in spite of large local effects. In the simplest case 



of two resonances, Eq.(31) would read 





2sin(A0 w )sin 


{A6(s 1 ,s 2 ))\A 2 Ai 
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2 + 2\A 2 A X 


cos(A0 w/ )cos(A5(si,s 2 )) 



A C p[Si,s 2 ) = — — — — — .,, A , W ., AS ., „ (5 



The asymmetry is driven by A8(si,s 2 ) due to interfering Breit-Wigner functions spread 
over the phase space. This is equivalent to having a different value of p for each bin. The 
distribution of A^p for the CP violating bins therefore depends strongly on the final state, 
on which resonances and with which relative phases it is built of. 

An important effect is the charge asymmetry induced by different production mech- 
anisms. This is not possible in pp collider, but it may occur in asymmetric collisions 
(like for LHCb data). The production asymmetry may as large as a 1-2% effect. Since 
it depends on the heavy meson momentum, it may vary across the Dalitz plot. In the 
following examples we assume that any eventual production asymmetry would lead only 
to a global effect, constant throughout the Dalitz plot. 
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4.1 Comment on CPT Constraints 



It is mentioned usually that CPT symmetry gives equality of masses and total widths of 
P and P. However it gives also equality of different classes of final states, where mixing 
happens; some general comments are given in Sect. 4.10 in (TS]. Up to isospin violation 
one has for example: 

F(B uAs -»■ 2vr, KK, 4vr, 2K2K, 6vr) = T(B uAs -»■ 2n, KK, 4vr, 2K2K, 6n) (59) 
T(D u>d>a -> 2tt, 4tt) = r(D U)d) . -+ 2tt, tftf, 4tt) (60) 

While mixing happens - and diagrams show it - we have little quantitive control over 
it. In a qualitative way one expects correlations like between CP asymmetries in D° —> 
K + K~ and D° — > 7r + 7r~ or in Bd — >■ K~7T + and Bd — > Kgit etc. As emphasized before 
that CP asymmetries with three-body final states will give us more information - and 
probably crucial one - about the underlying dynamics. Since SM produce sizable CP 
violation in b — > sqq with q = u,d,s one expects sizable CP asymmetries in Bd — > 
Ksp°/Ks<j/K + p~ /kit and higher resonances etc. with different signs compensate for 
CPT relation - but only qualitatively in practice. 

4.2 'Miranda Procedure' for Bd Three-Body Decays 

We give 'realistic' studies for Bd — > Ks7t + tt~ , where we have decent data and some infor- 
mation about the resonant structure [THl dZj. In all studies we consider time integrated, 
tagged samples. In each exercise two samples of B°,B° — > Ksir + ft~ were simulated inde- 
pendently using the same set of resonant amplitudes, namely Ksp, Ksfo(9S0), i£s/o(1370), 
K*(892)n and K$Xc- The samples are generated with CP violation seeded in three dif- 
ferent ways, as described above. 
A few remarks are in order: 

• Indirect CP violation has been very well measured in B d — > ipK$ with sin20!//3 = 
0.658 ±0.024; this observable enters in many transitions as an input quantity. How- 
ever for the time integrated Bd + Bd rates indirect CP asymmetry cannot contribute 
- unless there is a production asymmetry for Bd vs. Bd- 

• Direct CP violation can occur even in the time integrated Bd + Bd rates. 

• In the SM one has three quark-level processes, namely two tree-level b — > uus and 
b — > suu, where the second one is generated by QCD radiative corrections, and the 
loop Penguin b — > s + g's. They produce another dd and uu pair for the final state 
Kstt + it~ and a ss for Bd — > K$K + K~ . The Penguin operator b — » s + g's generates 
no weak phase; since it produces a AJ = transition, there is no appreciable 
relative strong phase from this contribution. On the other hand b — > uus and 
b suu represent a combination of AI = and AI = 1 amplitudes that in 
general will have different strong phases. As an example for ND: Charged Higgs 
exchanges would probably affect mostly b — > uus and b — > suu, introduce another 
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weak phase and different strong phases. Furthermore they should affect final states 
with pseudoscalar-scalar more than pseudoscalar- vector. The impact of ND in direct 
CP violation should be clearer in the former than the latter, since the latter 'suffer' 
from a larger 'background' from CKM. 

• Sizable contributions from final states K s p, K*(892)ir, i^*(1430)7r, K s fo(980) and 
i£s/o(1500) have been reported. No obvious contributions from Ks<J and/or kit 
have been found, but might be hidden under the 'no-resonance' listing without 30 
% of the rate of B d — y K s tt + it~. There are several theoretical arguments that such 
final states Kgo and kit should occur in an appreciable way. 

• While CKM dynamics has been found to produce at least the leading source of 
indirect CP violation in B d — B d oscillations, ND could still represent up to about 
20 % of it. While no clear deviation from CKM theory has been found in direct CP 
asymmetries in B d decays, ND could produce significant contributions. One expects 
that the weight of ND in direct CP asymmetries will change differently for classes 
of channels like pseudoscalar-vector vs. pseudoscalar-scalar. 

4.2.1 B d /B d ->■ K s n + n- - Constant CPV 

Direct CP violation has been found in B d — y K + n~ around 10 %. No sign has been found 
in C(B d — y K n°) = 0.00 ± 0.13. Yet one could find sizable impact with future data. 

The first and simplest study of the 'Miranda Procedure' refers to the case where one 
has one single source of direct CP violation acting on a given region of the Dalitz plot. 
The CP violation is seeded as a 10% excess of B° over B in the region s^ s7r +, 5^+^- < 7.5 
GeV 2 /c 4 . We have generated 300K B° and 330K B decays dividing the combined Dalitz 
plot into 256 bins of equal population. This excess of -B over B° events is equivalent to 
a global A C p of 4.76%. The distribution of the A^p across the Dalitz plot is shown in 
Fig.l. 

Having only one source of CP violation (constant p), the values of A^p for the bins in 
the region where CP violation was seeded are the same within statistical fluctuations. We 
therefore expect the A^p for the CP violating bins to be also distributed as a Gaussian. 
The distribution in Fig.2 is fitted by two Gaussian functions. The one representing the CP 
conserving bins has fixed mean (// = 0) and sigma (a = l/y/N), whereas the parameters 
defining the second Gaussian are free. 

The average value of Aqp in the region where CP violation was seeded is the mean of 
the second Gaussian, (13.64 ± 0.25)%. The normalization of each Gaussian is the number 
of bins that conserve/ violate CP. There are 167±13 bins conserving CP and 89±9 bins 
in which CP is violated. The number of events is the same for all bins, so we can obtain 
the global Acp from the ratio of CP violating to the total number of bins, and from the 
average value of Aqp, 

A CP = 712 < A c { p >= 4.98 ± 0.54% (61) 
n\ + n 2 
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Figure 1: Distribution of A^p 1 across the _B,j — > Kstt + tt~ Dalitz plot. In this example 
a single source of CP violation - constant excess of Bd over Bd restricted to the low 
K s7r + / 'tt + ti~ mass region - was simulated. Bins have different size in order to contain the 
same number of events. 

We not only recover the global Aqp but also access the average A^p and the fraction of 
events that violate CP and thus the localization of the source. 

This exercise clearly shows how the relatively large local effect is diluted when the CP 
violation strength is measured by the global Aqp (in this case, by the ratio of the area of 
the CP violation region and the total Dalitz plot area). 

4.2.2 BdjBd — > K s n + ii~ - Difference in Magnitudes 

In this example CP violation is seeded as a 10% difference in the magnitude of the resonant 
mode K*(892)n, a x , (892)7r = 0.9aK^(892)n. 

The total decay rates of B° and -B are proportional to the integral over the phase space 
of Ai and A4, respectively. In the present example, this means a global CP asymmetry of 
2.1%. Samples of 300K B° and 287K B decays were generated. The combined B° and 
B° Dalitz plot was divided into 1024 bins of equal population. 

The extra B° events are distributed in the bins along the if* (892) band, and the 
resulting A^p across the Dalitz plot is shown in Fig. 3; note that the values of A^p in the 
CP violating region are always positive. 

In Fig.4 the distribution of A^p for all bins is presented. The bins with no CP violation 
have equal number of B° and B° decays, within statistical fluctuations, resulting on a 
Gaussian distribution of A^p with jj, — and a = 1/y/N. The distribution of A^p for 
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Figure 2: Distribution of A^p for single, constant source of CP violation. The distribution 
was fitted to two Gaussians. The one centered at zero represents bins where CP is 
conserved, whereas the second Gaussian represents the CP violating bins. 

the CP violating bins is again parameterized by Gaussian, but this is used as an effective 
representation (p is no longer constant). As in the previous example, the distribution in 
Fig.4 is fitted to two Gaussian functions, one with fixed mean and sigma representing the 
CP conserving bins. 

We find 838 ± 46 bins with no CP violation and 186 ± 38 bins with average value 
of A^p 1 (11.1 ± 1.7)%. From these parameters we extract the global asymmetry, Aqp = 
2.1 ±0.1%. 

These exercises show that the observable A^p 1 carries the relevant information about 
local asymmetry. In both cases the CP violation is restricted to certain regions of the 
Dalitz plot, but leads to a global asymmetry. The local effects are much more intense 
than the phase space integrated ones. We showed that the later can be recovered in a 
consistent way. 

4.2.3 Bd/Bd Ksn + i^~ — Difference in Relative Phases 

We now discuss a more general case, where the CP violation occur via a difference between 
B and B in relative phase of a given set of resonances. This is a much more difficult and 
subtle situation, which depends strongly on the final state characteristics: (i) Which 
resonances are present, (ii) What are their relative phases, (iii) Is there an sizable 
contribution from scalars [IS] . 

Two independent samples were generated using the same set of resonances as in the 
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Figure 3: Distribution of A^p across the Bd — > Ksk + 7[~ Dalitz plot for the case of CP 
violation seeded as a difference in the magnitude of the pK$ mode between and B^. 
The excess of B^ over Bj events is concentrated along the p band. Note that the values 
of A]jp for these bins are always positive. 
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Figure 4: Distribution of A^p for the case of CP violation through a difference in the 
magnitude of the pKs mode. The Gaussian in red is an empirical representation of A^p 
for the CP violating bins. 



previous examples. A 60° phase difference in the pK$ mode was introduced between 
Bd and Bd- The combined Dalitz plot was then divided into 1024 bins. The seeded 
phase difference is large, causing local asymmetries that can be as large as 80%, shown 
in Fig.5. The global asymmetry, however, is small: 1.0%. Due to the phase variation of 
the Breit-Wigner curve, the CP asymmetry change sign along the p band. In this case 
the integration over the phase space - necessary to compute the total rates - cancels out 
most of the effect of CP violation. This cancellation is clearly seen in Fig. 6, which has an 
enlarged view of the Dalitz plot region where CP violation occur. 

The interference between pK$ and the other resonant modes, which is governed by 
the combined strong phases of the Breit-Wigner curve, A<5(si, s 2 ), causes each bin to act 
as an independent source of CP violation. When one has repeated the same experiment 
many times, the values of A^p for each bin would be distributed with a mean and sigma 
), respectively, each bin having its own value of p. The distribution 
of A^p for all bins will have two components, as in the other examples. The distribution 
from the CP violating bins would no longer be a Gaussian, but some function that is 
particular to each specific final state. 

In general there would be as many bins with positive and negative ^4cp±> so the P ro ~ 
cedure adopted in the previous examples would underestimate the measurement of the 
average asymmetry in this case. As before, we can fit the distribution to a Gaussian for 
the CP conserving bins - \i — and a = 1/yfN, but with unknown area - plus one 



given by Eqs.d54|55 
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Figure 5: Distribution of A^p across the Bd — > Ks^ + ^ Dalitz plot for the third example. 
When the CP violation is seeded as a relative phase difference, the values of A^p change 
sign. 



function for the CP violating bins. Having defined the Gaussian CP conserving bins, this 
can then be subtracted off,and two numbers could be computed: the average value of 
^4cp-i- for the regions where the asymmetry is either positive or negative. 

We illustrate this procedure in Fig. 7, where the A^p distribution for the CP violating 
bins was empirically fit to two Gaussians. The fit yields 568 ± 83 bins in the CP conserving 
Gaussian, and 466 ± 41 bins in which CP is violated. We then compute the weighted 
average value of ^cpV ^ or ^ ne ne S a ti ve an d for the positive part of the distribution in 
Fig. 7. This yields 

< A%£_ >= -(14 ±2)%, (62) 

and 

< >= (16 ± 2)%. (63) 

In order to test this procedure, we go to the limit of very high statistics. Since the 
width of the distribution of CP conserving bins is a = 1/yN, when iV is very large, 
the Gaussian gets very narrow, in practice restricted to the central bin of Fig. 7. We can 
then compute the weighted average in the negative and positive regions separately with 
a simple counting procedure, discarding the central bin. The average values of ^cpt 
obtained are 

< A^_ >= -15% (64) 

and 

< A^ + >= 22.2% (65) 
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Figure 6: Enlarged view of the CP violating bins of Fig.5. 




Figure 7: Distribution of A^p for the third example. In addition to the Gaussian represent- 
ing the CP conserving bins (green curve), two other Gaussians were used to empirically 
represent the AqP distribution of the CP violating bins (red and blue curves). 
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in good agreement with the fitting procedure used in the more realistic scenario. 

5 B s Three-Body Decays 

At present, the experimental situation is very different for B s transitions even beyond the 
fact that B s — B s oscillations are very fast. 

• Within SM indirect CP violation is small - i.e. sin2/3 s ~ 0.03 — 0.05. Finding it sig- 
nificantly larger is a clear manifestation of ND. There is some evidence that indirect 
CP violation is larger than predicted by CKM; studies of B s — y ipcj) in LHC data 
should clarify this issue and allow sin2/3 s as an input for searching manifestations 
of ND. 

• If one indeed finds that CKM theory does not produce the leading source of indirect 
CP violation, there is a good chance that ND generates a larger contribution also 
to direct CP violations. 

• Measuring y s more accurately will help in cross checking finding CP asymmetries in 
the sum of time integrated B s and B s rates. 

• There are no data for B s — y K s tt + it~ or B s — y K S K + K~ . One expects the Dalitz 
plots for these B s transitions very different for these Bd transitions. 

• The tree diagram b — y uud and the Cabibbo disfavoured penguin one-loop diagram 
b — y d + g's to generate direct CP violation in both CKM and ND. Again final states 
like K s cr should show clearly manifestations of the impact of ND. 

6 Outlook 

Present data from Belle, BaBar, CDF and LHCb and future ones from LHCb, Super-Belle 
and Super-BaBar have reached the status to probe the possible impact from ND in B u ^, s 
and D u ^, s with accuracy and correlations. Analyzing non-leptonic three-body final states 
there needs significantly more experimental efforts through 'Miranda Procedure' - but it 
will be awarded with more lessons about the underlying dynamics and deep insights into 
its 'shape'. 

The Miranda Procedure I is a good way to show whether or not there is CP asymmetry 
in three-body decays of D and B mesons. It can also tell us where in the Dalitz plot CP 
violation occurs and give hints of the kind of operators that are involved. A further 
development of this technique, presented here, is a necessary step towards a quantitative 
output. One should keep in mind, however, the crucial difference between two- and 
three-body decays: while in former case CP asymmetries are observed in total decay 
rates, in the latter case there are several options for CP violation manifestations. CP 
asymmetries through phase difference - 'favorited' by model builders - are intrinsically 
complicated because each bin acts as an independent source of CPV. Moreover, strong 
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interactions governing phases across the Dalitz plot are still out of control quantitively. 
Accurate data on three-body final states will help efforts from theorists working on HEP 
and Hadrodynamics/MEP. 

One should not forget about constraints from CPT symmetry, but those are not of 
quantitative value on a practical level; it should tell us to think about other channels in a 
qualitative level; theoretical inputs help here. 

The 'Miranda Procedure' helps greatly to 'localize' CP asymmetries and find evidence 
for the impact of ND and its 'shape'. It does not mean that theoretical inputs are not 
needed, but to focus on them. It should enhance interests from theorists working in HEP 
and HP/MEP. 

Applying 'Second Generation of Miranda Procedure' is now at the 'starting line' - the 
'race' will proceed over many longer 'distances' with simulations and - most importantly 
- with real data: 

• Time integrated and non-flavour tagged rates for B UtS /D U} d, s decays; 

• Flavour tagged ones for B uAs /D uAs ; 

• partially time integrated ones; 

• t — y u[K tt/K2tt/?>K] decays. 

One needs no more hardware - 'only' thinking and working. 
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